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A few quasi-exactly solvable models are studied within the quantum Hamilton-Jacobi formalism. 
By assuming a simple singularity structure of the quantum momentum function, we show that the 
exact quantization condition leads to the condition for quasi-exact solvability. 



I. INTRODUCTION 



This paper reports an investigation of quasi-exactly solvable (QES) models within the quantum Hamilton-Jacobi 
(QHJ) formalism. One of the earliest investigations on QES models was done by V. Singh et.al [1] . The QES models 
have been studied extensively and for a general review and references we refer the reader to the book by Ushveridze 
[2]. In the recent years, a lot of interesting work has been done on QES periodic potentials. We refer the reader to 
references [3-13]. A complete list of canonical forms for the QES models has been obtained in [3,4]. The QES models 
have a characteristic property that, when the potential parameters satisfy a specified condition, analytic expressions 
for a few energy levels and their corresponding eigenfunctions can be obtained exactly. This condition between the 
' parameters of the potential will be referred to as quasi-exactly solvability condition. 

A well known example of QES model is the sextic oscillator in one dimension, corresponding to the potential 
V(x) = ax 2 + fix 4 + 7X 6 . The condition of quasi-exactly solvability is found to be "^[§^ — a ] — 3 + 2n, where n is 
a non-negative integer which is related to the number of levels for which exact energy eigenfunctions and eigenvalues 
\Q , can be computed. A number of other QES models have been constructed and studied within algebraic and group 
CO ■ theoretical approach. Interesting connections have been established between the sextic oscillator and second order 
linear differential equations within a new approach to second order linear differential equations [13]. 

In this paper we study the sextic oscillator and a few other QES models in one dimension. In each of these cases 
the condition of quasi-exactly solvability is derived within the QHJ approach. The next section contains an overview 
of QHJ formalism and how it is used in the present paper to study the QES models. In the following four sections 
we investigate the sextic oscillator, sextic oscillator with a centrifugal barrier, a hyperbolic potential and a circular 
*£h potential. The last section contains our conclusions about exact solvability and quasi-exactly solvability. 
P-l. 
^— > ' 

c ■ 

=5 

The QHJ formalism was initiated by Leacock and Padgett [14,15] and was successfully applied to several exactly 
solvable models (ESM) in one dimension in quantum mechanics by Bhalla et.al [16-18]. Our discussions in this paper 
•'■h will be limited to the QES problems in one dimension only. The Schrodinger equation is 

On substituting tp = exp(iS/h), in (1) we obtain the following equation for S 

f) 2 -«0=M^-^)]. (2) 
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II. QUANTUM HAMILTON-JACOBI FORMALISM 



dx J \ dx 2 

We define p(x) = which satisfies the Riccati equation 



p 2 (x) - ih-^-p(x) = 2m[E - V{x)\. (3) 
dx 
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In the limit 7i — > 0, (2) becomes the classical Hamilton- Jacobi equation and p(x) — > p c i, the classical momentum which 
is 



p cl = y/2m[E-V(x)\. (4) 

Therefore (2) and (3) will be referred to as QHJ equations, p(x) will be called the quantum momentum function 
(QMF) and S will be the quantum action. In terms of the eigenfunction ip of the energy, the QMF is given by 

p W - . (5) 

An important step in the QHJ formalism is to regard i as a complex variable and to extend the definition of p(x) to 
the complex plane. From (5) it is obvious that the zeros of the wavefunction correspond to the poles of the QMF. 
It is known that the bound state solutions of (1) corresponding to the n th level has n real zeros correspondingly the 
QMF has n poles on the real line. From (3) it can be seen that if a; is a point at which V{x) is analytic and p has 
a pole, the pole must be of first order and the residue at that pole will be — ih. Therefore the integral of the QMF 
taken along the contour C, which encloses these poles, will have the value nh. Thus we get 

— (f p{x)dx = nh. (6) 
27r Jc 

and this is an exact quantization condition in one dimension [14,15]. This exact quantization condition (6) has been 
used to obtain bound state energy levels without solving for the cigenfunctions for several ESM [16-18]. For this 
purpose one needs to know the location of singularities and their corresponding residues of the QMF in the complex 
plane. The QMF has the two kinds of singularities, fixed and the moving singularities. The fixed singularities corre- 
spond to the singularities of the potential and will be present in every solution of the Rlccati equation, their location 
being independent of the initial conditions. The position of the moving singularities depends on the initial conditions. 
It is known that for the solutions of a Riccati equation only poles can appear as moving singularities. Therefore if 
the potential is meromorphic, the solutions will also be meromorphic. Coming back to the QMF corresponding to the 
solutions of the Schrodinger equation, in addition to the poles, corresponding to the n real zeros of the n th excited 
state, in general there may be other moving poles. The knowledge about these poles is needed to apply QHJ method. 
The residue at any fixed pole can be computed from (3). This, being a quadratic equation, will lead to two solutions 
and therefore a boundary condition on the QMF is needed to pick up the physical solutions. Leacock and Padgett 
proposed that one should make use of the condition, 

\imp(x)^p c i. (7) 

h— >o 

From (4) it should be noted that the classical turning points become the branch points of p c i . It has to be emphasized 
here that, when a; is a complex variable, p c i is a double valued function. A full definition of p c i in the complex plane 
as a function of a complex variable x, requires us to select a particular branch of p c i which is assumed to have a 
branch cut in the classical region and a positive value just below the branch cut. Locating the singular points of the 
QMF and imposing the boundary conditions discussed above are two crucial but usually difficult steps in arriving at 
the correct solutions. In the earlier studies of ESM, it was easy to guess the singularity structure of the QMF and it 
was found that, there were no moving poles away from the real axis. For the QES models studied in this paper, it is 
very difficult to find the location of the singularities of the QMF. In order to be able to make some progress we make 
a simplifying assumption about the moving poles and the nature of the singularity at infinity. The details of moving 
and possible fixed poles will be given in each case separately as and when we discuss the model. We now state our 
main assumption, common to all models studied, i.e, the point at infinity is an isolated point and that it is a pole of 
finite order and not an essential singularity. Under these assumptions we show that, imposing the exact quantization 
condition leads to the condition of QES in each case. 



III. SEXTIC OSCILLATOR 



The potential for the sextic oscillator is 

V(x) = ax 2 +f3x 4 + ~/x 6 , 7 > 0. (8) 

From now onwards we set h = 2m = 1. For the n th exited state, QMF has n moving poles on the real axis and we 
assume that there are no other moving poles. In order to use the quantization condition (6) we need to evaluate the 
integral J(E) 
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J(E) = -L / p{x)dx (9) 



along the contour C. As a first step, it is convenient to deform the contour C to a large circular contour T, with the 
center at the origin and big enough to enclose all the singularities of the QMF in the finite complex plane. This is 
possible because the point at infinity is assumed to be an isolated singularity and the value of the integral remains 
unchanged. Hence 

J(E) = ±- £p(x)dx. (10) 

In order to evaluate the integral along the large circle T, we apply an inversion mapping x — > y = —. Under this 
mapping the circular contour T gets mapped to another circular contour 7 in the y plane. The only singular point 
inside this contour is the point y = 0, which is the image of the point at infinity in the x-plane. Thus we obtain from 
(10) 

J[E) = h i m 7 dy - 

To evaluate the integral in (11) a Laurent expansion of p(y) is made. 

3 7 oo 
n=l y n=l 

Substituting this expansion in (11) and integrating term by term, we get 

J{E) = iai. (13) 

It only remains to compute the coefficient a\ of the Laurent expansion given in (12). To do this we start from the 
QHJ equation 



P 2 {y) + iy 2 ^p{y) = E-^-^-^. 

dy y 2 y 4 y fi 



Substituting the Laurent expansion and equating the coefficients of different powers of y in both sides of the equation 
(14), we successively get the following equations. 

b 2 3 = - 7 . (15) 
6 2 = 0. (16) 
26i6 3 - -/?■ (17) 

b\ + 6 3 (2ai - 3i) = -a. (18) 

It is important to know that, we would get two solutions for 61 corresponding to the two solutions of 63 = ±1^/7. 
This happens due to the fact that the QHJ is quadratic in the QMF. Thus one needs a boundary condition to pick 
the correct solution. We propose to use the square integrability of the wavefunction instead of the original boundary 
condition, explained in the introduction, which was proposed by Leacock and Padgett. This is because the original 
boundary condition is difficult to implement in the present case due to the presence of six branch points in the p c \. 
In order to find the restrictions coming from the square integrability, we compute the wavefunction 

4>(x) = cxp ip(x)dx^j (19) 

for large x as follows. The most important term in the Laurent expansion (12) for small y w 0, corresponding to large 
x, is 
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m « h fs (20) 

and the wavefunction for large x becomes 

(63X 4 \ 
i - J dx. (21) 

Out of the two solutions, 63 = ±1^/7, -0 is square intcgrable only for 63 = Using this value of 63 and calculating 
J(E) from (13), the quantization condition gives 

In order to compare the result in (22) with the known condition [2], we write 

7 = a 2 , /3 = 2ab. (23) 
Thus we get a = b 2 — a(3 + 2n) which agrees with known result. 

IV. SEXTIC OSCILLATOR WITH THE CENTRIFUGAL BARRIER 

This kind of model was first studied by A.P.Houtot [19]. The potential we discuss here is 

V(x) = 4(5 - - + [b 2 - 4a(5 + l - + M)\x 2 + 2abx i + a 2 x 6 (24) 

and the range of S is taken to be 45 > 3. We note that the potential goes to 00 as x — > and x — > 00. Hence the 
classical turning points will be on the positive real axis. The physical motion in both classical and quantum situations 
will be confined to positive real axis only. As, has been discussed earlier there will be n moving poles in the classical 
region which are enclosed by the contour C in the quantization condition stated in (6). When we extend the definition 
of x, to take all the complex values, we expect n additional moving poles on the negative real axis. These come from 
the symmetry of the potential under the transformation x — ► —x. Just as in the quantization condition (6), we have 

<p p(x)dx = n (25) 

where C\ is the contour which encloses the n additional moving poles on the negative real axis. The QHJ equation 
for this potential is 

p 2 (x) - i^-p(x) = E - A(S - ]){S - h\ - [b 2 - Aa{S +\ + M)\x 2 - 2abx 4 - a 2 x 6 . (26) 

We observe that in addition to the 2n moving poles, x = is a fixed pole. In this case we assume that there are no 
other singularities in the finite complex plane. As in the case of the sextic oscillator, in order to evaluate J(E) in (9), 
we deform the contour C to a large circular contour T, which will enclose the 2n moving poles and the fixed pole at 
x = 0. Therefore we have, 

&> p(x)dx = &> p(x)dx + <x> p(x)dx + (p p(x)dx (27) 
Jr J~f Jc Jd 

where 70 is the contour enclosing only the fixed pole at x = 0. Expanding p{x) in the Laurent series, in the powers 
of x 

, 00 

P(x) = ^ + Y,anx n (28) 

n=0 

we get 

-!- * p(x)dx = ib x . (29) 
2tt 7 70 
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We fix b\ by substituting the Laurent series in (26). Thus we obtain, 

6i = ^[4S-3], &! = -i[4S-l]. (30) 
The choice for the value of b\ consistent with square integrability in the specified range is 

6 1 = -i[4S-l]. (31) 

Thus (27) becomes 

— / p{x)dx = ibi + 2n. (32) 

The integral along T is computed by changing the variable from x — > y = ^ and proceeding in the same way as was 
done for the sextic oscillator. From (31) we then obtain 

M = n. (33) 

which is the quasi-exact solvability condition for the sextic oscillator with the centrifugal barrier. 



where, 



V. CIRCULAR POTENTIAL 

The potential for this model is 

A B 

V(x) = — 5— + = C sin 2 x + D sin 4 x (34) 

sin x cos 2 x 

A = 4(Si - -\),B = A(S 2 \)(S 2 I) (35) 

C = q\ + 4«i(Si + S 2 + M), £> = q\ (36) 
and the ranges of 5i and S 2 are 2S\ > 1 and 25*2 > 1. The QHJ equation is 

p 2 (x) - i 4-p(x) =E ^ 4— + C sm2 x — D sin 4 x (37) 

ax sin x cos 2 x 

and the quantization condition is 

— (b p(x)dx = n (38) 

Define sin 2 x = t, the QHJ equation becomes 

p 2 (t)-2 l ^(T^T)j t p(t) = E-j--^- t +Ct-Dt 2 (39) 



where p(t) = _p(x). Defining o by p = yt(l — t)q we obtain the QHJ equation in t variable as follows, 

2 9 d i(l-2t)q _ E A B C Dt 

q ~ Jt q+ t(l-t) ~ t(l-t) ~ t 2 (l-t) ~ t{l-tf + ~ T { ' 

The quantization condition (38) becomes 

<?(*) 



We observe that QMF has fixed poles at t = 1,0 and proceeding in a way similar to §4, we obtain the quasi-exact 
solvability condition as 

M = n. (42) 
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VI. HYPERBOLIC POTENTIAL 



The QHJ equation for this potential is 



where 



p 2 (x) - i -^-p(x) = E H ^ ^— + C cosh 2 x-D cosh 4 x (43) 

ax cosh x sinh a; 



A = 4(^ - -\),B = A(S 2 \){S 2 \) (44) 



C = q\ + 4«i (Si + S 2 + M), D = q\ (45) 

and the ranges of S\ and 5 2 arc 2S\ > 1 and 2S 2 > 1. A change of variable from a; — > £ = coshx changes 
p(x) — > p(t) = p(cosh _1 a;), then defining 

p(t) = Vt 2 -l9(t). (46) 

The QHJ equation becomes 

2 _.d_ £_ B | Ct 2 gfo 4 

9 ^ 9 t 2 -l t 2 - 1 t 2 (t 2 - 1) (t 2 -l) 2 t 2 -l t 2 -l 1 j 

The quantization condition is 

±-Iq(t)dt = n (48) 

We expect 2n moving poles on the entire real line as in Sec. IV, and fixed poles at t = 0, ±1. These poles together 
with a pole of finite order at infinity lead to the condition 

M = n. (49) 

which agrees with known condition of QES of this model. 



VII. CONCLUSIONS 



In the limit h — > QMF of the sextic oscillator goes to p c i which has six branch points. Therefore, in general 
we expect a complicated singularity structure for the sextic and also for other models studied in this paper. In the 
analysis presented in the previous sections it was assumed that there are no moving poles off the real axis. A closer 
look at the derivation of the condition of quasi-exact solvability shows that the assumption" no moving poles of the 
QMF off the real axis " can be replaced by a weaker assumption, " QMF has a finite number of moving poles in the 
complex plane " without altering any of the results. In fact it can be seen from explicit solutions in [2,4] that the 
algebraic eigen functions do have complex zeros. For all the models that have been studied here, the algebraic part 
of the spectrum and eigenfunctions are well known and in all the cases QMF has a pole at infinity. Thus we conclude 
that, for the class of potentials studied here, QES models are the only models for which QMF has a pole at infinity and 
a finite number of moving poles in the complex plane. For each model the known algebraic eigenfunctions correspond 
to QMF having singularity structure as postulated. For completeness it must be mentioned that the assumption of 
finite number of poles is not independent of the assumption that the point at infinity is a pole. For a large class 
of potentials which are analytic everywhere, except for isolated singularities, the moving singularities of solutions of 
QHJ can only be poles. An infinite number of such poles will, therefore, have an accumulation point at infinity and 
making it z = oo a non- isolated essential singular point. 

The parameter n that appears in the exact quantization condition is related to the number of moving poles in QMF 
and has different roles to play for ESM and QES models. In ESM each value of n corresponds to an energy level 
and an cigenfunction with n real zeros. In case of the QES models, n appears as a parameter in the expression for 
the potential and picks out a particular QES model within a family of potentials ; varying n gives rise to a different 
potential within the family. Recalling that n also determines the number of moving poles of the QMF it appears 
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reasonable to expect that all the algebraic eigenfunctions for a given QES models (fixed n) will have the same number 
of the complex zeros determined by n. An explicit check reveals that this expectation is true for the sextic oscillator. 
A preliminary study reveals that, the poles of the QMF for QES periodic potentials have a richer structure. A detailed 
study of the location of the poles of the QMF in different QES models will be reported elsewhere. 
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